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Topologically ordered quantum systems have robust physical properties, such as quasiparticle 
statistics and ground-state degeneracy, which do not depend on the microscopic details of the Hamil¬ 
tonian. We consider topological phase transitions under a deformation such as an effective string 
tension on a Z 3 topological state. This is studied in terms of the gauge-symmetry preserved quan¬ 
tum state renormalization group, first proposed by He, Moradi and Wen [Phys. Rev. B 90 , 205114 
(2014)]. In this approach modular matrices S and T can be obtained and used as order parameters 
to characterize the topological properties of the phase and determine phase transitions. From a 
mapping to a classical 2D Potts model on the square lattice, the critical string tension, at which 
the transition to a topologically trivial phase takes place, can be obtained analytically and agrees 
with the numerically determined value. Such a transition can be generalized to a Z n topological 
model under a string tension and determined in the same way. With different deformations, the 
Z3 topological phase can also be driven to a critical phase which contains, in the large deformation 
limit, the wavefunction analogous to the Rokhsar-Kivelson point in the quantum dimer model in 
one case and the fully packed loop model in another case. 


I. INTRODUCTION 

Conventional phases and their transitions in 
condensed-matter systems can be understood by 
Landau symmetry breaking [T|. However, new phases 
have emerged over the past few decades such as integer 
quantum hall [2j [3] and fractional quantum hall [4] 
effects have evaded the usual Landau-Ginzburg-Wilson 
paradigm. More recently, microscopic spin models 
exhibiting topological orders were also constructed, 
such as the toric code, quantum double models and 
string-nets [SHE]. These new quantum phases are called 
topological phases and they cannot be characterized by a 
local order parameter. Instead, they are found to possess 
the so-called topological order 0 uni characterized by 
properties such as the ground-state degeneracy, nontriv¬ 
ial quasiparticle statistics EHEl, and more recently 
nonzero topological entanglement entropy mm- 

Recent progress on the general scheme of classifying 
these topological orders is also vitalized by using the 
notation of entanglement [181122] . Specifically, intrinsic 
topological orders have patterns of long-range entangle¬ 
ment that cannot be achieved via local unitary trans¬ 
formation from a trivial product state. The local uni¬ 
tary transformations can only remove short-range entan¬ 
glement, i.e., quantum correlations between neighboring 
sites. Due to its robustness against the local operations, 
the long-range entanglement gives rise to topological or¬ 
der of the nontrivial ground state. Long-range entangle¬ 
ment is manifested in the so-called topological entangle¬ 
ment entropy nu [ 13123 H 2 H 1 El] and other related entan¬ 
glement quantities [291132] . 

However, in general, it is still numerically challenging 
to access large system sizes to extrapolate accurately the 
topological entanglement to determine whether a topo¬ 
logical phase exists. Even if one manages to obtain 
nonzero topological entanglement entropy, there are pos¬ 
sibly different topological phases that possess the same 


value of the entanglement. 

Quasiparticle statistics, related to the generally non- 
abelian geometric phases, represented by the modular S 
and T [33] matrices, provide a more refined characteriza¬ 
tion of a topological order. S and T are representations 
of 90°-rotation and the Dehn twist, respectively, on a 
two-dimensional torus. These topological quantities do 
not depend on the microscopic details of the Hamilto¬ 
nian, and can be used as order parameters, albeit non¬ 
local. It turns out that one can, given a complete basis 
states in the degenerate ground space, exploit entangle¬ 
ment, with respect to a bipartite cut, to deduce the set 
of minimally entangled states (MESs). These states form 
a special complete set of bases [I5j [34l[38j in the ground 
space, and from them the modular matrices S and T can 
actually be deduced. Methods such as Quantum Monte 
Carlo and the exact diagonalization have been used in 
such computations [33 El 0o]. 

The recent development of matrix product states 
(MPS) [411143] and the 2D tensor product states 
(TPS) [441149] . which are the generalization of 
the density-matrix renormalization group (DMRG) 
method ISBEIl has given rise to alternative approaches 
to obtain ground states and thus MESs. These tensor- 
network or tensor-product states also were constructed 
as variational wavefunctions for optimizing Hamiltoni¬ 
ans Hasans G 21021 , but they have also been shown to 
exactly represent a large class of topological states, in¬ 
cluding both non-chiral [541157] and chiral [58, 5£] topo¬ 
logical order. The DMRG itself has also been used to 
obtain the MES systematically (aUEHUHOl- In particular, 
via techniques introduced in Refs. [35, 36? 39, 60, 61], the 
MESs corresponding to different quasiparticle excitations 
can be obtained and the modular matrices can be deter¬ 
mined from MESs. We shall follow the approach pro¬ 
posed by He, Moradi and Wen EH> who showed that the 
gauge-symmetry structure of TPS can give us informa¬ 
tion about topological order. In particular, applying the 


2 


gauge-symmetry preserved tensor renormalization group 
(GSPTRG) to Z 2 topological orders under deformation 
of the wavefunction via a string tension, He et al. ob¬ 
tained the modular matrices as a topological order pa¬ 
rameter to characterize the topological phases and their 
phase transition from the Z 2 topologically ordered phase 
to a trivial phase (that is adiabatically connected to a 
product state). Here, we employ the same approach and 
apply it to Ztv topological models under deformations 
and characterize the topological phase with the modular 
matrices and locate the phase transitions to topologically 
trivial phases. For one particular type of deformation, we 
can map the wavefunction norm square to the partition 
function of the classical AT-state Potts model and obtain 
the analytic critical string tensions. The numerical re¬ 
sults via GSPTRG agree very well with those via the 
mapping. The unexpected result we obtain is that un¬ 
der different deformations, the Z3 topologically ordered 
phase can also be driven to a critical phase, in addition 
to the trivial product-state phase. The phase diagrams 
of the Z 3 topological phase under different deformations 
considered in this paper are summarized in Fig. [14| 

The paper is organized as follows. In Sec. [TT] we re¬ 
view the notion of topological order and gauge symme¬ 
try preserved tensor renormalization group which can be 


used to identify intrinsic topological orders; In Sec. Ill 


we discuss how the deformation can be applied to the 
Z n topologically ordered model via the tensor product 
states and how to use GSPTRG to compute the modular 
matrices for characterizing the topological order. There, 
we also show a useful mapping form the Z^v model under 
a string tension to the two-dimensional classical TV-state 
Potts model, from which the critical point between the 
topologically ordered phase and a trivial phase can be 
obtained analytically and compared with numerics. In 
Sec. |IV[ we evaluate the modular matrices and correla¬ 
tion function for the Z3 model and discuss the critical 
phases. We conclude in Sec. [V] 


II. MODULAR MATRICES AND GAUGE 
SYMMETRY PRESERVED TENSOR 
RENORMALIZATION GROUP 


Specifically, to obtain the modular matrices, we need 
to first determine all the degenerate ground states 
{\^Pa)}a=i °f system. The S and T matrices are de¬ 
fined as follows mmm-- 

^ a \S\^ b )=e- asV+oi - 1 ^S ab 
^ a \f\i> b ) = e~ aTV+0 ^T ab , ( 1 ) 

where S and T are the transformations of the 90° rotation 
and the Dehn twist respectively on a torus with lattice 
size V , as and ax are non-universal constants, and S a b 
and are elements of the modular matrices. The in¬ 
formation of quasiparticles statistics and their fusion rule 
are encoded in the S and T matrices [33 J [631465] . 

In particular, from the gauge structure of tensor prod¬ 
uct states at the fixed point (TPS) [ 20 ) [ 66 ][ 68 ] , the de¬ 
generate ground states can be obtained by inserting the 
gauge transformation to TPS, i.e., operating on the bond 
or virtual degrees of freedom by the appropriate gauge 
transformation. The degenerate ground states can be la¬ 
beled as 1 2 p(g,h)) with gauge transformations (g,h) ap¬ 
plied on the internal indices along two directions. 

Next, we describe the real-space renormalization group 
approach and how the modular matrices can be obtained. 


B. Gauge symmetry preserved tensor 
renormalization group (GSPTRG) 

An symmetry preserved tensor renormalization group 
procedure exist for quantum states based on the tensor 
product representation 

H) = Y tTv(A^A s \..A^...)\ Sl s 2 ...s m ...), ( 2 ) 

Sl,S 2 ,...S m ... 

where A® ^ . is a local tensor with physical index s 

and internal indices af3 7 etc. tTr denotes tensor contrac¬ 
tion of all the connected inner indices according to the 
underlying lattice structure. The norm of TPS is given 

by 


(V#> = tTr(T 1 T 2 T 3 ...T m ...), (3) 


A. Modular matrices 

One of the exotic features of the topological order is 
the nontrivial quasiparticle statistics, which can be ob¬ 
tained by modular transformations on degenerate ground 
states on the torus, giving rise to the modular S and T 
matrices. This approach of characterizing topological or¬ 
ders has become quite fruitful recently fl5l l36l fBTl [62] . 
The modular matrices, or S and T matrices, are gener¬ 
ated respectively by the 90° rotation and the Dehn twist 
on a set of degenerate ground states on the torus. The 
elements of S matrix express the mutual statistics of the 
quasiparticles, whereas the T matrix express the twisting 
a quasiparticle wavefunction along an axes by 360°. 


where the local double tensor TP can be formed by merg¬ 
ing two layers tensors A and A* with the physical indices 
contracted, 

T = Z( A «,/3.7,a-) X (4) 

However, it is generally computionally hard to calculate 
exactly the tensor trace (tTr) or the contraction of the 
whole tensor network in two and higher dimensions This 
imposes the hurdles of an exponentially hard calculation. 
Several approximation schemes have been proposed as so¬ 
lutions in this context, such as iPEPS [S3 algorithm, the 
corner transfer matrix method (CTMRG) [69! , and ten¬ 
sor renormalization approach |48[, all of which tackle this 
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problem essentially by scaling the computational com¬ 
plexity down to the polynomial level for calculating the 
tensor trace. 

We express tensor renormalization group (TRG) ap¬ 
proach which is akin to the real space renormalization in 
the way that at each step, the RG is structured by merg¬ 
ing sites (by contracting respective tensors) and truncat¬ 
ing the bond dimension according to the relevance of the 
eigenvalues in the Schmidt decomposition of the old ten¬ 
sors. Actually, by doing several steps of TRG, the double 
tensor will flow to a fixed point. The information of the 
ground states can be extracted from the fixed point ten¬ 
sor. 

A many body wave function and Hamiltonian may 
be invariant under certain transformations that corre¬ 
spond to symmetries. The symmetry group divides the 
Hilbert space of system into symmetry sectors that can 
be labeled by the quantum number or conserved par¬ 
ticle number. Symmetries can be used to improve the 
numerical methods, such as used in the exact diago- 
nalization method and density matrix renormalization 
group E2BIZD, as well as matrix product states and ten¬ 
sor product states [6T1. 72, 73] . The TRG method has it¬ 
self emerged as an alternative approach for dealing with 
quantum spin systems. However, the local decomposition 
usually breaks the symmetry sectors if the symmetry is 
not strictly enforced. Then, it would flow to the wrong 
fixed point. In this section, we will express the gauge 
symmetry preserved tensor renormalization group (GSP- 
TRG) procedure proposed by He, Moradi and Wen [6l] 
for TPS and demonstrate its effectiveness in identifying 
the topological order from the fixed point tensors. GSP- 
TRG differs from TRG only when we decompose the ten¬ 
sor. The matrix would has blocks corresponding to the 
symmetry sectors. Once we keep the block structure, the 
symmetry will be preserved. 



FIG. 1. Schematic procedure of the GSPTRG 

In Ref. [51] , the authors showed how to preserve the 
Z 2 gauge symmetry to determine the phase diagram of 


Z 2 and double-semion model with string tension. The 
Z 2 topological model can be realized in toric code model. 
The ground state of toric code model is an equal weight 
superposition of all closed string loops and it can be rep¬ 
resented by tensor product state with virtual dimension 
X = 2. When putting the Hamiltonian on the torus, it 
has four-fold degeneracy ground states which can be ob¬ 
tained by inserting the string operators, such as Pauli 
matrix Z, to the virtual bond degrees of freedom of one 
particular ground-state wavefunction. 

Below we shall consider the case of Z^v topologically 
ordered phases for which the Hamiltonian is the gener¬ 
alized toric code model HU [ 75 , so the Zjv topological 
order has Z^v gauge symmetry. To implement GSP¬ 
TRG, first we form a double tensor T on each site 

T = T, s ( A a4 3,7,5) X i A a',6')* aS shoWn in FigR 1 )' 
The double tensor T will has Z^ x Z n gauge symmetry. 
Due to such a gauge symmetry, the elements of double 
tensor are non-zero only when a + /3-\- / y + S = 0 (mod 
N ) and a' + /3' + 7 ' + 5' = 0 (mod N). 

We then view the tensor T as a matrix M a/ 3 a f^ :1 sY8' = 
T 0 : 07 / 3 / 3 ', 77 ',< 5 <v- The first step of the coarse graining is to 
decompose a rank-four tensor M (e.g., the double tensor 
T) into two rank-three tensors. We do it in two different 
ways on black and white tensors (see Fig. [l] (3)). Due to 
a such gauge symmetry, the tensor 

N 

M-otficx 1 /3' , 7 < 57 ' 8' = ( 5 ) 

p,q =i 

would be block diagonalized by the quantum number. 
For example, the each block m p , q obey the rule 

a + /3 = p mod N 

7 + S = N — p mod N 

a' + /3' = q mod N 

j' + S' = N — q mod N , ( 6 ) 

where p,q = 0,1, 2 , 3,..., N — 1 as shown in Fig. 1 (b)- 
Then, singular-value decomposition (SVD) is performed 
in each block. As mentioned above, the tensor contrac¬ 
tion is an exponentially hard calculation. Here, a cut¬ 
off x c might be necessary on the dimension of double 
tensor to keep the computation efficient. When making 
the truncation, we need to preserve the symmetry struc¬ 
ture of the tensor by keeping the blocks together. Such 
symmetry considerations apply similarly to all symmetry 
groups. 

After the decomposition, the lattice structure is 
changed. The second step is to form a new rank-four ten¬ 
sor denoted by T' as shown in Fig. [l] (5). To do this, we 
combine the resultant four tensors that meet at s square 
to form a new tensor as shown in Fig. [T] (4). After doing 
several steps GSPTRG, the double tensor will gradually 
flow to a fixed point tensor that preserves the gauge sym¬ 
metry. 

The modular matrices can be evaluated and monitored 
during the process of the RG steps by performing three 
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FIG. 2. (a) The double tensor structure, (b) The M is block 
diagonalized by the quantum number. 


steps, (i) Inserting the gauge transformations g,h,g',h' 
into the internal indices a,/3, a',/?', respectively of the 
fixed point double tensor as shown in Fig. [ 2 ] (a) to deter¬ 
mine 

(ii) Performing the rotation and the Dehn twist oper¬ 
ators on the ground state wave function. The operators 
performing on the physical indices can be achieved or 
replaced by appropriate gauge operations to internal in¬ 
dices, 

('ip(g',h , )\T\'ip(g,h)) = (^(g 1 , h')\ip(g,gh)) 
(ip(g',h')\S\ip{g,h)) = (ip(g' ,h')\^{h,g~ 1 )). 

(iii) Finally, tracing all internal indices of the fixed 
point tensor. 

For the Z 2 topological phase in Ref. m, the Z 2 gauge 
symmetry is generated by cr z acting on each internal in¬ 
dices. (We note that for different forms of ground-state 
wavefunctions, the gauge operator might be a x .) The 


same rule holds in the Zjv topological phase that gener¬ 
alizes from this case. The gauge symmetry can be gen¬ 
erated by the TV x TV operator Z at which all elements 
are zero except diagonal term Zk k = exp 2?r2 ^~ 1 ) • k = 
1, 2,3, ...N. 

Instead of viewing the gauge operators g , h applied to 
internal indices, we can also understand the degenerate 
ground state \^p(g,h)) in terms of closed string opera¬ 
tors < 7 , h defined in the physical dimension. The S and T 
transformation can be written in the basis of the degen¬ 
erate ground states |^(g, ft)) generated by appropriate 
closed string operators g,h on the torus along the ver¬ 
tical and horizontal directions, respectively. The string 
operators for the particular Z^r phase (i.e. Z^v toric code) 
can be chosen as Z q = ( Z) q 0 (Z^) q 0 {Z) q 0 {Z^) q ..., 
q = 0,1, 2, ..., TV — 1 and Z° = 1. The |^(X,X)> = \^) 
is our reference ground state represented by the tensor 
product state in Eq. [ 8 ] 

Our goal is to obtain the modular S and 
T matrices. We can obtain the TV 2 x TV 2 
^-matrix generated by 90°-rotation and it is 
given as S a , b = (iP{Z a \ Z a2 ))\S^(Z b \ Z b2 )) = 

(‘ip(Z ai ,Z a2 )\‘ip(Z b2 ,Z~ bl ))] where ai,a 2 ,bi,b 2 = 
0,1, 2, ...TV — 1. Both matrix indices a and b range from 
1 to TV 2 where a = a\ TV + <22 + 1 and b = bi TV + 6 2 + 1. 
Similarly, to apply the Dehn twist to all degener¬ 
acy ground states we also can get the TV 2 x TV 2 
T-matrix: T a , 6 = {^{Z a \ Z a ^)\T\^{Z b \ Z b ^)) = 

(ijj(Z ai , Z a2 \ijj(Z hl , Z hl + b2 )y For example, the 4-by-4 
T-matrix for the Z 2 toric code is given as follows: 


<^(x,x)i^(x,x)) (ip(i,i)mi,z)) (^(i,i)mz,z)) (ip(i,i)^(z,i)) \ 
U(i,z)\i>(i,i)) U(i,z)\^(i,z)) U>(i,z)\^z,z)) U(i,z)\^(z,i)) 
(^(Z,1M(1,1)) (iP(z,i)mi,z)) (ip(z,imz,z)) 

\(ii>(Z,Z)\il>(I,l)) (ii>(Z,Z)\il>(I,Z)) (ti>(Z,Z)\il>(Z,Z)) (ip(Z, Z)\ip(Z,X)) ) 


III. THE RESULTS 


A. The quantum Zjv phase 


Let us begin by describing the construction of Z^v 
wavefunctions. Their Hamiltonian is generalized from 
toric code model heej m - The tensor product state 
(TPS) on the square lattice motivated by the Z^r topo¬ 
logically ordered phase is characterized by the rank-4 
tensor, P a ,f 3,^,6 with four internal indices running over 
0,1, 2 , ...N — 1 on vertex and the rank-3 tensor G s a p with 
one physical index s running over the TV possible spin 
states 0,1, ....(TV — 1) on the link as shown in Fig. [ 6 ] (a). 


The wave function is then given by 

H) = F tTr (®t>-P G Si ) I Si, S 2 , (8) 

{*’} 

where v labels vertices and l links. Specifically, 

p _ / 1, if a + /3 + 7 + S = 0 mod N, , . 
“’ /3 ’ 7 ’ <5 - I 0j otherwise, W 

and 

GU = 1, 2 = 0,1,2,3,..., TV-1, 

others= 0 . ( 10 ) 

The rank-3 tensor G behaves like a projector which 
essentially sets the internal index equal to the physical 
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FIG. 3. The trace of modular matrices (a)S', (b) T, and the property (c) X 2 /X 1 as functions of parameter g display a phase 
transition at critical point g c of Z 3 model 


index. In m, they studied the problem of the stability 
of a tensor network state under physical perturbations 
to the local tensor. In view of this we can consider a 
deformation, Q = qi\i)(i\ and 0 < g* < 1 which 

apply to the physical indices, \^(q)} = Q 0 Q 0 ... 0 
Ql'ip). At qi = 1, (i = 0,1, 2,..., N — 1), this is exactly Z^v 
topologically ordered phases. At go = 1, q^ = 0, (i = 
1 , 2 , 3 , ...,7V — 1 ), the tensor represent a product state of 
all 0. At some critical point in parameters g^, the phase 
transition will occur. Here, we first consider go = 1 and 


Qi = g 2 ,{i = 1 , 2 , 3 ,..., TV - 1 ). 

The Ztv phase has a TV 2 -fold ground-state degeneracy 
on a torus, which corresponds to N 2 different types of 
quasiparticle excitations. Therefore, the corresponding 
modular matrices will be of size N 2 x TV 2 . For simplic¬ 
ity of the calculation we associate every vertex with four 
matrices as shown in Fig. [6] and from the double tensor. 
The norm of wave function represented by the double ten¬ 
sor can then be represented as standard tensor product 
form as Eq. ( 13 ), where the double tensor f ^ has 
eight inner inaices a/, /?', 7', S ', <a, /?, 7, S = 0 , 1 ,..., N — 1 . 
From the TPS with deformation, we can find the phase 
transition point of the Z3 model as shown in Fig. [U b y 
using GSPTRG. The fixed point tensor structure might 
be complicated but it is always possible to identify them. 
We calculate S and T matrices and a basis independent 
quantity given by the ratio X2/X1 [ 2 F, where X\ and 
X2 as shown in Fig. [ 4 ] are defined as follows, 


Xl = ( £ A U<°,P x t4a',/3',a',/3')*F 

s,Oi,(3,a' ,/3' 

X 2 = (^a,/3, 7 ,/3 X ^- 7 ,<5,a,<5) X 

s,s ', a, f3 , 7 ,5, a' ,/3' , 7 ' ,8' 

({Ai',p>,y,p>r x (11) 
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FIG. 4. The quantity X 2 /X 1 obtained by taking the ratio 
of the contraction value of the double tensor in two different 
ways. X 2 /X 1 is invariant under gauge transformation, such 
as unitary operators U and V. It can be used to distinguish 
different fixed-point tensors. 


0.7776 < g < 1 . 0 , the tensor belongs to the Z 3 topologi¬ 
cally ordered phase, since we obtain nontrivial S and T 
matrices as follows: 


S = 


/l 

0 

0 

0 

0 

0 

0 

0 

\o 


0 0 
0 0 
0 0 
1 0 
0 0 
0 0 
0 1 
0 0 
0 0 


0 0 0 
0 0 0 
1 0 0 
0 0 0 
0 0 0 
0 1 0 
0 0 0 
0 0 0 
0 0 1 


0 0 
1 0 
0 0 
0 0 
0 1 
0 0 
0 0 
0 0 
0 0 


°\ 

0 

0 

0 

0 

0 

0 

1 

0 / 


( 12 ) 


T = 


/l 

0 

0 

0 

0 

0 

0 

0 

V0 


0 0 0 
1 0 0 
0 1 0 
0 0 0 
0 0 1 
0 0 0 
0 0 0 
0 0 0 
0 0 0 


0 0 0 
0 0 0 
0 0 0 
0 1 0 
0 0 0 
1 0 0 
0 0 0 
0 0 0 
0 0 1 


0 0 \ 
0 0 
0 0 
0 0 
0 0 
0 0 
1 0 
0 1 
00 / 


( 13 ) 


We find that when 0 < g < 0 . 7776 , all components 
of S and T matrices are 1, and X2/X1 = 1.0, and this 
shows that the ground state is in the trivial phase. When 


and X2/X1 = 0 . 33333 . The S and T matrices obtained 
above give us the modular transformations. Note that 
the matrices are not in the canonical form (where, e.g., 
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the T is diagonal), but there is a procedure to make the 
T-matrix diagonal and at the same time make S in the 
canonical form [64] . Then diagonalized T-matrix gives 
the self-statistics of quasiparticles, and S matrix the mu¬ 
tual statistics. 

In the GSPTRG algorithm, we preserve the Z 3 gauge 
symmetry of the tenser. As performing more steps of 
RG, the nonlocal order parameters show sharper changes 
around g c = 0.7776 in Fig. [3j and the crossings of differ¬ 
ent RG curves signal the transition point separating the 
trivial phase and topological phase. 

Similar behaviors of the Z4 and Z5 model under GSP¬ 
TRG is also found in Fig. [5] For Z4 model, when 
0.7597 < g < 1.0, the tensor belongs the Z4 topologi¬ 
cally ordered phase, since we obtain nontrivial S and T 
matrices as follows: 


/! 


S = 


1 


(14) 


/I 


T = 


1 / 


(15) 


and X 2 /Xi = 0.25 as shown in Fig. 5] (a). For Z 5 model, 
when 0.7450 < g < 1.0, the tensor belongs the Z5 topo¬ 
logically ordered phase, since we obtain nontrivial 25 x 25 
S and T matrices and X 2 /X 1 = 0.2 as shown in Fig. p] 
0 >). 

In the following, we shall find the mapping from the 
norm square of the Zat wavefunctions and the TV-state 



(a) 



FIG. 5. The trace of modular matrices S as functions of 
parameter g display a phase transition at critical point g c of 
(a)Z 4 (b) Z 5 model 


Potts model. First, by applying a deformation Q = 
Eto 1 to tensor G® a new tensor A can be ob¬ 

tained as shown in Fig. ©(b): 

A if = £ Qs,s'G s ' p X Q:, s „g£. (16) 

s,s' ,s" 


There are two nonzero components in A 

Aoio = 1; Ki = Qi (* = 1, 2 ,N — 1 ). (17) 

Second, we form the double tensor T as shown in Fig. [6] 
(e) from tensors P = P 0 P* on the vertex and A on the 
link as shown in Fig. ©( d )> 


r w'i' ,j' ik' ,1' 


E 


Pa,g,j,6 x ,P',^',8' 


a, ( 3 , 7 , 8, a', ( 3 ' ,7', 8' 



m =0 


where rim is the number of virtual indices in state ”ra”. 
Then the double tensor with nonzero components of Z 2 
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FIG. 6. (a) The tensor product state representation of Zjv 
topologically ordered phase, (b) Apply the deformation to 
form a double tensor A on the link and the tensor on the 
vertex, P — P <g> P *, (c) The double tensor represented by A 
and P, (d)Decompose tensor A = and reconstruct the 

tensor (e) The new double tensor. 


B. The classical iV-state Potts Model 

The Hamiltonian of the TV-state Potts is given by 

#=E(- J W- ( 22 ) 

<i,3> 


with the sum running over the nearest neighbor pairs 
<i,j > over all lattice sites. The degree of freedom s*, Sj 
on the site is on values in {0,1, — 1}. The S Si , Sj is the 

Kronecker delta, which equals one whenever Si = Sj and 
zero otherwise. The model is ferromagnetic when J > 0 
and antiferromagnetic if J < 0. A basic question is where 
the phase transition point is. Baxter has determined the 
exact free energy for the square-lattice Potts model and 
determined the critical point, such as ferromagnetic crit¬ 
ical point is e^ J = 1 + y/q and antiferromagnetic critical 
point is e^ J = — 1 + ^4 — q (76] . 

For example, the three-state Potts model with zero 
external field is N = 3 and then the spins are usually 
taken to be {0,1,2}. By tuning the temperature, the 
phase transition will occur, for example, for J > 0 the 
transition is first order if N > 5 and is continuous if 
N < 4. 

The partition function of TV-state Potts Model is given 

by: 


model for Q = | 0 )( 0 | + g 2 |l)(l| are given by 


npOOOO _ -i 

0000 — 1 

npOOll _ npOHO _ nr. 1100 

Mon — Miio — Mioo 


nrlOOl _ npOlOl 

1001 ~ MlOl 


_ tpIOIO _ 4 

— Moio — 9 


nnllll 

Mm 


= 9 


( 20 ) 


The double tensor with nonzero components of Z3 model 
for Q = | 0 )( 0 | + g 2 |l)(l| + g 2 |2)(2| are given by 


o o 
o o 
o o 
o o 

= 1 





o o 

o o 

II 

o o 
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The Z 2 model above is mathematically equivalent to 
the two-dimensional classical Ising model where the tran¬ 
sition point is known to great accuracy [Tfi . The numer¬ 
ical results from Ref. ED show that the Z 2 toric code 
and double-semion models have the same critical point 
under the string tension. This is understood by recogniz¬ 
ing that the phases in the double-semion wavefunction of 
Ref. [ 61 ] under the string tension cancel in the mapping 
to the Potts model and its has the same norm square as 
the one from deforming the toric code. Therefore, they 
have the same critical point g c = 0.802243 from Eq. ( 31 ) 
below. Next we derive the mapping. 


z=e-? H =Y J n ePJSsi - sj =e n ( 23 ) 

<i,3> <i,3> 


Here, A is named the transfer Matrix on the link and can 
be represented by a N x N matrix 


A (i,j) 


(e 0J 1 1 ... 1 \ 

1 e 0J 1 ... 1 

y 1 1 1 ... e 0J ) 


( 24 ) 


Applying the Hadamard matrix H, we form a diagonal 
matrix A' as shown in Fig. [ 7 ] (a). 


A' =H x A x H ] 

fe 0J + q- 1 0 0 ... 0 \ 

0 e 0J -l 0 ... 0 

\ 0 0 0 ... e 0J -1 J 


( 25 ) 


For example, the Hadamard matrix of three-state Potts 
model can be given by 


2 i /! 1 1 \ 

H = = - 7 = 1 <*> b ( 26 ) 

a,( 3=0 V 3 y ^ CJ 2 UJ ) 


where u = e* 277 / 3 . 



















FIG. 7. The tensor representation of partition function, (b) 
Reconstruct the tensor (c) The new tensor configuration. 


Then, we define a trivial rank-4 tensor R a aaa = 1 with 
index a running over 0,1,2, ...,7V — 1 on each vertex. To 
apply the Hadamard matrix H to the tensor R , 

Ri,j,k,l = Rg,a,a,a x jH+ k H a ,l 

a 

f 1/q, if<a + /3 + 7 + £ = 0 mod N 
1 0 , otherwise. 

Combining the tensor R' and A' forms a new tensor T. 
Then the double tensor with nonzero components of clas¬ 
sical Ising model without magnetic field are 


(27) 

(28) 


separateing the ordered (ferromagnetic) and disordered 
(paramagnetic) phases. The form of partition function of 
the TV-state Potts model represented by the tensor net¬ 
work is equal to the norm square of the Zw wave function, 
with deformation Q = |0)(0| g 2 \i) (i\, represented 

by the tensor network, such as Eq. ([2Q|, Eq. (29) and 


Eq. (36), Eq. (30). The double tensor of the norm of Zw 


wave function is just a two copies of the partition func¬ 
tion of Potts model. We then have the relation between 
the parameter g and the parameter /?J, 


9 = 


VeP J - 1 
\fePJ + N - l 2 


1/8 


(31) 


From the relation and the transition point of Potts 
model, we can then obtain the phase transition point g c 
for Ztv model : 



Numerics 

From mapping 

Numerics 

Potts model 

N 

9c 

9c{PJ) 

1 /" 


2 

0.8021 ( Xc = 24) 

0.802243 

1.010 

1 

3 

0.7776 (xc = 24 ) 

0.777817 

1.201 

6/5 

4 

0.7597 (xc = 32 ) 

0.759835 

1.501 

3/2 

5 

0.7450 (xc = 30 ) 

0.745582 


first order 


Toooo = ^(2 cosh (/?J)) 2 

Toon = Tono = Tnoo = -(2 cosh(/ 3 J ))(2 sinh(/ 3 J)) 

Tiooi = Toioi = Tioio = -(2 cosh(/ 3 J ))(2 sinh(/ 3 J)) 

Tim = I( 2 sinh(/ 3 J)) 2 . ( 29 ) 

The nonzero components of the three-state Potts 
model are given by 

Toooo = ^{Ve /3J + 2) 4 

Tom = Tion = Tnoi = Tmo = ~{V^ J + 2 )(\/e /3J — 1 ) 

T0222 = T2022 = T2202 = T2220 = + 2 )(\/e /3jr — l ) 3 

T0012 = T0120 = T1200 = T 2 ooi = 

T0021 = T0210 = T 2 ioo = T1002 = 

T0102 = T1020 = T0201 = T0201 = -(a/ e /3J + 2) 2 (V e /3J — 1) 
T1122 — T1221 = T2211 = T2112 = 

T1212 = T2121 = I(\/e^ - l) 4 . ( 30 ) 

C. Duality 

The two-dimensional classical ferromagnetic Potts 
models have phase transitions located at e^ cJ ) = 1 + A /q, 


From this table, the transition points from the GSPTRG 
are quite close to exact mapping results. 

For the classical TV-state Potts Model, if TV < 4, this 
model describes a continuous phase transition with scalar 
order parameter. The critical exponents of these transi¬ 
tion are universal values and characterize the singular 
properties of physical quantities. In particular, the cor¬ 
relation function of the classical spin system is, 

D{fi - fj) = (, S(n)S(fj )) - (S(fl))(S(fj)) ) (32) 

where the brackets mean statistical average over all con¬ 
figurations. The correlation length £ is defined in terms of 
correlation function D(fi — fj ) ^ e \ r i- r j\R^ where \ fl — fj\ 
is the distance between two spins. In the asymptotic limit 
of large \ f\ — fj |, the correlation function decays to zero 
exponentially. On the other hand, the correlation length 
diverges at the critical point £ ^ \T — T c | -zy , where v is 
an example of a critical exponent. In the above Table, 
we list the exponents for the classical Potts Model. 

As shown above the Zw model with one parameter 
g can be mapped to the classical TV-state Potts Model. 
We fit critical exponent v for the order parameter i.e. 
Tr(T),Tr(S), or X 2 /X 1 under the renormalization flow. 
Our results as shown in Fig. [ 8 ] are \jv = 1.01 for the 
Z 2 model, \/v = 1.201 for the Z 3 model and \/v = 
1.501 for the Z4 model. For the five-state Potts model, 
the transition is a weak first-order |77j. Near critical 
point, the correlation length is very large but finite. It 
is difficult to distinguish the critical exponent of the Z 5 
model numerically. It is worth mentioning that since 
after each step of the renormalization, the number of sites 
is reduced by half, this is to say that, after performing 
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FIG. 8 . The trace of S matrix as a function (g — gc ) x L u on (a)Z 2 (b)Z 3 (c) Z 4 models. 


two steps of renormalization, the distance between two 
new neighbor sites will be increased twofold. We can 
perform a rescaling, and see the data collapse. In Fig.[ 8 j 
L is length scale defined by the number of RG step n rg , 
L = 2 n ^/ 2 . 


IV. THE CRITICAL PHASE 


Here we focus on the Z 3 topologically ordered phase 
under different forms of deformation. As mentioned 
above, the Z3 topologically ordered phase can be repre¬ 
sented by rank-four tensor P (see Eq. ©) and rank-three 
tenors G (see Eq. © ). We have already discussed a spe¬ 
cial projector Q a = lx | 0 )( 0 |+g 2 x |l)(l|+g 2 x | 2 )( 2 | with 
0 < g < 1 in earlier parts, and here only briefly review it. 
At some critical point in g, the state must go through a 
phase transition. Moreover, this model is mathematically 
equivalent to two-dimensional three-state Potts model. 
We can obtain g c = 0.777817 from exactly mapping. By 
using GSPTRG, we can obtain the transition point to 
be between 0.7776 and 0.7777 and find our results to be 
with 0 . 02 % accuracy. 

We shall discuss the more deformation operator Q = 
q 0 x | 0 )( 0 | +q 1 x |1)(1| + q 2 x | 2 )( 2 | with qo,qi,q 2 > 0 in 
detail. By applying it to the Z3 toric code wavefunction, 
we obtain the double tensor 


Qo° Qi 1 Q 2 2 > if i+j + k + l = 0 mod 3 
0 , otherwise, 

(33) 


where i,j,k,l = 0 , 1,2 and no, ni, and n 2 means the 
number of the inner indices in state ” 0 ”, ” 1 ”, and ” 2 ” 
respectively. 


A. For q 0 = 1 \qi= 0; q 2 g 2 


We first study the case qo = l;qi = 0;q 2 = g 2 , i.e., 
the wavefunction is constructed from an effective two- 
dimensional Hilbert space spanned by |0) and |2). At 
g = 0 , the tensor represents a product of all states being 
0. Therefore, the phase diagram near g = 0 is a region 
of the trivial phase that is adiabatically connected to a 
product state. At g > 0 , the nonzero components are 


rrp0222 _ 

0222 ~ 


6 rrp2022 _ 6. np2202 _ 6. rrp2220 _ i 

iJ j ^ 2022 — y 5 ^ 2202 — i/ 5 ^ 2220 — £/ 


np0000 

0000 


= 1 . 


(34) 


While g 1, the tensor form is mathematically equiva¬ 
lent the quantum dimer model at Rokhsar-Kivelson (RK) 
point [ 8 Ql [ 8 l| [84] where the degenerate ground state at¬ 
tains exactly zero energy in the form of an equal weight 
superposition of all possible configurations in a given 
winding parity sector on square lattice. This point is 
a critical phase with algebraically decaying correlation 
function. It is noteworthy that the RK point is a criti¬ 
cal point separating two gapped phases of the quantum 
dimer model on the square lattice. In the large g limit, 
we can regard string -2 as vacuum state and string -0 as 
dimer. Then the tensor form will obey the hard-core con¬ 
straint where each vertex is connected to one dimer only 
and has the same weight. 

As a result, as g goes from 0 to 2 , a phase transition 
must occur. Some phases between gapped trivial states 
and critical phases could exist. Can a topologically or¬ 
dered phase exist? Let us examine this middle parameter 
region in detail. 

We use the gauge symmetry preserved tensor renormal¬ 
ization group to flow our wave function to fixed point and 
obtain the modular matrices by inserting gauge transfor¬ 
mation and invariant quantity X 2 /Xi easily with sharp 
quantum phases transition point. Again, this is a advan¬ 
tage to use GSPTRG to characterize the topologically 
ordered phase numerically with TPS ansatz. Our nu¬ 
merical results is given in Fig. [9j 
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FIG. 9. The trace of modular matrices (a)S', (b) T, and the property (c) X 2 /X 1 as functions of parameter g display a phase 
transition at critical point g c of Z3 model with projector Qb = lx |0){0| + 0 x |1)(1| + g 2 x |2)(2|. 


We see that when 0 < g < 0.944, all components 
of 9 x 9 S and T matrices are 1 , and X 2 /Xi = 1 . 0 . 
We can say the ground state in the trivial (product) 
phase, since at g = 0 this wave function is product 
state. For 0.944 < g < 1.238, the tensor belongs the 
Z 3 topologically ordered phase, since we obtain nontriv¬ 
ial S and T matrices (same as Eq. ( 12 ) and (13)) and 
X 2 /X 1 = 0.3333. For 1.238 < g < 2.0, when renor¬ 
malization step is larger enough, we can obtain the triv¬ 
ial fixed point. The GSPTRG can distinguish different 
topologically ordered phases and a topologically ordered 
phase from topologically trivial phases. But it cannot 
distinguish different topologically trivial phases. From 
the equivalence to the RK point for large g , the system 
is probably gapless for g large; from the trivial product 
state at g = 0 , the system is probably also a trivial prod¬ 
uct phase for small g. 

The GSPTRG can be used to detect the topological or¬ 
der well, but it cannot use to classify the trivial phases. 
It must look more carefully in trivial phases. We need 
sufficient evidence to show what trivial phase is. The rel¬ 
evant quantity of our interest is the correlation function 
which appears useful in characterizing the product state 
and critical phase. The connected correlation function is 
expressed as, 



D{r) = (So(f i )So(f i )) - (So(rj)){S 0 (rj)), (35) 

where r = |r* — fj\ and S'o(g) count 1 if "0" state is 
present on the link fi in a given configuration otherwise 
0. In the asymptotic limit of large r, the correlation 
function will converge to zero. 

When 0 < g < 0.944, we obtain topological entan¬ 
glement entropy 7 = 0. In the limit g = 0, this is a 
product state of all 0. There is good evidence to show 
that as 0 < g < 0.944 the state in this region is the trivial 
state. When g > 1.238, based on the observation of ten¬ 
sor representation on g 1 discussed earlier, we would 
claim that this phase might be a critical phase. From the 


FIG. 10. The correlation function under the deformation 
Q h = 1 x |0)(0| + 0 X |1)(1| + g 2 x 12)<21 with (a) g = 1.02 
(b) g m 2.0. The correlation function under the deformation 
Qc = 1 x |0)(0| + 1 x | 1 )( 1 | +g 2 X 12)(21 with (c) g = 1.3 (d) 
9 = 2 . 0 . 


GSPTRG results, the values of tr(S), tr(T) : and X 2 /Xi 
approach to 9.0, 9.0, and 1.0 respectively. In the tensor 
renormalization approach, the above quantities are cal¬ 
culated on the square lattice with size 128 x 128. For the 
critical phase, the correlation function is algebraic of the 
form D(r) ~ r~ b for large r as displayed in Fig. 10 (b) 
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g g g 

FIG. 11 . The trace of modular matrices (a)S', (b) T, and the property (c) X 2 /X 1 as functions of parameter g display a phase 
transition at critical point g c of Z3 model with projector Q c = lx |0)(0| + 1 x 1 1) (1 1 + g 2 x |2)(2|. 


obtained using the tensor renormalization group (TRG) 
and the mean-field approximated second renormalization 
group (SRG) [79] . However, for TRG, the convergence 
becomes very slow and unstable in the critical phase [78]. 
Even though the TRG is not efficient, the results with 
large critical bond dimension Xc = 64 already show that 
the correlation function structure is different exponen¬ 
tially decaying as shown in Fig. 10 (a). On the other 
hand, the mean-field approach of SRG can improve the 
accuracy of the results as shown in Fig. [lO] (b) with bond 
dimension Xc = 64. When 0.994 < g < 1.238, it is worth 
noting that the Z 3 phase exists, even if one of physical 
term |1) is turned off. It is an unsettled question. 


B. For q 0 = 1; q ± = 1; q 2 g 2 


Here we consider the deformation parametrized by 
q 0 = 1 ,qi = 1, q 2 = g 2 in the projector Q c = q 0 x |0)(0| + 
qi x 11)(11 + g 2 x |2)(2|. The double tensor with nonzero 
components of are given by 
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At g = 1, this is exactly the Z 3 phase and the corre¬ 
sponding state has topological order. Recall the earlier 
case, the same may be said, while g 1, the tensor form 
is also mathematically equivalent to the quantum dimer 
model at Rokhsar-Kivelson (RK) point on the square lat¬ 
tice, the wavefunction at which is the equal weight su¬ 
perposition of all dimer configurations. At some critical 
point in g , the phase transition will occur. 


With the GSPTRG procedure, we flow the initial dou¬ 
ble tensor for arbitrary g to symmetry preserved fixed 
point tensor and calculate the modular matrices and the 
X 2 /X 1 . The results are shown in Fig. 


11 


As the 


ber of GSPTRG steps increase, the transition from Z 3 
phase to a critical phase approaches to a step function 
at g c = 1.76. From the correlation function, we also see 
that it is algebraically decaying as shown in Fig. [To| (d). 


C. For q 0 = 1; q\ = g 2 ; <?2 = g 2 


In the earlier section, we only consider that g goes 
from 0 to 1. However, in the large g limit, the effective 
terms of Eg. [36| are T{™ = Tji{ = = T*™ = 


¥{212 = Tiili = 9 8 - We can regard string-2 as vacuum 
state and string-1 as dimer state. It will obey the rule of 
fully packed loop model where each vertex is connected to 
two dimers. Just as the RK point of the quantum dimer 
model on the square lattice, the wavefunction at the large 
g limit corresponds to the equal weight superposition of 
all fully packed loop. The RK-like point of fully packed 
loop model is also a critical liquid state with algebraically 
decaying correlation functions [82] . 

The transition from Z 3 phase to the critical phase oc¬ 
curs at g c = 7.3825 from GSPTRG. We also calculate 
the correlation function by using SRG, and see the alge¬ 
braically decaying beharior as shown in Fig. [12] 

The GSPTRG can be used to classify the topologically 
ordered phases with different modular matrices. The 
scheme is robust for topological order, because it is a 
gapped state and has gapped entanglement spectrum or 
singular value (SV) spectrum. When performing the RG 
transformation, as long as the cutoff in the bond dimen¬ 
sion Xc is large the wave function will flow to a fixed 
point. However, the scheme cannot be used to distin¬ 
guish different non-topological (trivial) phases. As we 
have observed, gapped and gapless trivial phases under 
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FIG. 12. The correlation function under the deformation 
Qa — lx |0)(0| + g 2 x 11)<11 + g 2 x |2)(2| with g = 20 .0. 


the GSPTRG scheme flow to fixed points that exhibit the 
same tr(S), tr(T), and X 2 /X 1 . However, this approach 
to fixed point is slower for the gapless phases as shown 
in Fig. [9] 

In Fig. [l3j we plot the spectra of singular values af¬ 
ter several steps of GSPTRG transformation. For g = 0.9 
which is a topological phase, the there-fold degeneracy of 
SV spectra can be found and the pattern of degeneracy 
is robust under RG flow. The single largest SV separated 
by gap can be found in the gapped trivial phase. How¬ 
ever, due to the finite cutoff y c , the RG flow for the gap¬ 
less phase cannot be followed accurately and indefinitely. 
The continuous spectrum of SV eventually breaks into 
chunks separated by a gap after sufficiently large num¬ 
ber of RG step; see Fig. 13 (b). Increasing the bond di¬ 


mension y c can slow the breakdown. An alternative, and 
perhaps the most elegant approach is proposed by Even- 
bly and Vidal [83] , who proposed a coarse-graining trans¬ 
formation called tenser network renormalization (TNR) 
that can explicitly recovers the scale invariance and flow 
the wavefunction to a fixed point for critical points. How¬ 
ever, implementing the TNR scheme is beyond the scope 
of the present work. 


D. Summary of phase diagram 

The results of the section suggest the generic Z 3 model 
with projector Q = q 0 x J0)(0| + qi x |1)(1| + q 2 x | 2 )( 2 | 
shown in Fig.[l4] For projector Q a , by varying parameter 
g , the phase transition from trivial state to topologically 
ordered phase will occur at g c \ = 0.7776. The second 
phase transition from trivial state to critical phase occur 
at g c 2 = 7.3825. For projector Q&, by varying parameter 
g , the first phase transition from trivial state to topo¬ 
logically ordered phase will occur at g c \ = 0.944. The 
second critical point from topologically ordered phase to 
critical phase is at g c 2 = 1.238. The interesting thing is 
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FIG. 13. Spectra of singular values of g=0.7, g=0.9 and 
g=15.0 for cutoff Xc = 24 after (a) 6 and (b) 20 applications 
of GSPTRG, under the deformation Q a = 1 x |0)(0| + g 2 x 
axil + p 2 x |2)(2|. The index i is to number the singular 
values A-s. 


that even we turn off one parameter, the Z 3 phase can be 
found with two body system. For projector Q c , by vary¬ 
ing parameter g , the phase transition from topologically 
ordered phase to critical phase will occur at g c \ = 1.76. 
In general, richer phase diagrams may be obtained by 
considering more general parameter q 0 , q \, and q^ or 7L^ 
model. 


E. For N > 3 case 

In general, for other Z^v phases, we always can find 
the phase transition with deformation Q = (ft N) I • 

From tensor structure, sometimes we can deduce the pos¬ 
sible phase. For example, we consider the Z 4 phase with 
deformation Q = | 0 )( 0 | + | 1 )( 1 | + | 2 )( 2 | + g 2 x |3)(3|. As 
g 1 , the tensor from represent a product state of all 
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= 1 x |0)(0| + g 1 2 x 11)(1| + fir 2 x |2)(2| 
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FIG. 14. Generic phase diagrams, (a) the deformations 
Qa = 1 x |0>(0| + s 2 x |l)(l|+5 2 X |2)(2|, (b) Q b = lx |0)(0| + 
Ox |l)(l|+p 2 x |2)(2|, (c) Q c = lx |0)(0| + lx |l)(l|+(? 2 x |2}(2|. 


3. This is because the T 3?3;3?3 was allowed. However, for 
the Z 3 phase, the 72,2,2,2 was not allowed. Thus, it is 
impossible to find the product state for Z 3 phase with 

Q = |0)(0| + 11)(l| + g 2 x |2)(2|. 


V. CONCLUIONS 


We have employed the gauge symmetry protected 
tensor renormalization group (GSPTRG) method intro¬ 
duced by He, Moradi and Wen Ref. |6T] to study Z^v 
topological order under deformation. It is important to 
know the underlying gauge symmetry operators in order 
for this method to work. Due to the removal of irrelevant 
short-range entanglement, the fixed-point wavefunction 
contains primarily the long-range entanglements that can 
be used to identify the topological order. From the fixed- 


point form of the tensor representing the ground-state 
wavefunction, the modular matrices S and T that repre¬ 
sent the mutual and self-statistics, respectively, of quasi¬ 
particles can be obtained. These can be used as order 
parameters to detect phase transitions. We applied the 
string tension to deform the Zjv wavefunction, similar to 
that in the Z 2 toric code with deformation. The GSP¬ 
TRG approach accurately determined the phase transi¬ 
tion between the nontrivial topologically ordered phase 
and the trivial phase, the result of which matches very 
well with the mapping to TV-state Potts model. The RG 
process is a coarse-graining process and can be associated 
with the change of a length scale. From this perspective 
we were able to collapse the data from modular matrices 
with suitable scaling near the transitions and determined 
the critical exponent of the correlation length and the re¬ 
sult agrees well with the mapping to the Potts model. In 
particular, we also investigated different deformations on 
the Z 3 model and we found that the topologically or¬ 
dered phase can be driven to a critical phase. Moreover, 
there exists a finite region of parameters such that the Z 3 
phase is composed of local two-level systems, i.e., qubits. 

There has been tremendous progress on both theoreti¬ 
cal and experimental advancement of the search of exotic 
phases with topological order. A future extension and 
application of the GSPTRG method would be to first 
use a gauge symmetry preserved approach to find the 
ground-state wavefunction of a Hamiltonian (potentially 
containing topological order) and use the RG to flow the 
wavefunction to a fixed point and obtain the modular 
matrices. 
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